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ABSTRACT 

We present the first high spatial resolution near-infrared direct and polarimetric ob- 
servations of Parsamian 21, obtained with the VLT/NACO instrument. We comple- 
mented these measurements with archival infrared observations, such as HST/WFPC2 
imaging, HST/NICMOS polarimetry, Spitzcr IRAC and MIPS photometry, Spitzer 
IRS spectroscopy as well as ISO photometry. Our main conclusions are the follow- 
ing: (1) we argue that Parsamian 21 is probably an FU Orionis-type object; (2) 
Parsamian 21 is not associated with any rich cluster of young stars; (3) our mea- 
surements reveal a circumstellar envelope, a polar cavity and an edge-on disc; the disc 
seems to be geometrically flat and extends from approximately 48 to 360 AU from the 
star; (4) the SED can be reproduced with a simple model of a circumstellar disc and 
an envelope; (5) withi n the framework of an ev olutionary sequence of FUors proposed 
bv I Green et al.l ()2006f ) and lQuanz et al.l ()2007f ). Parsamian 21 can be classified as an 
intermediate-aged object. 

Key words: stars: circumstellar matter - stars: pre-main sequence - stars: individual: 
Parsamian 21 - infrared: stars - techniques: polarimetric. 



1 INTRODUCTION 

The most important process of low-mass star formation is 
the accretion of circumstellar material onto the young star. 
According to current models, the mass accumulation is time- 
dependent with alternating episodes of high and low accre- 
tion rates (|Vorobvov fc Basul l2006l : iBolev et al.ll2006l ). The 
high phase may correspond to FU Orionis objects (FUors), a 
unique class of low-mass pre-main sequence stars that have 
undergone a m ajor outburst in optical light of 4 mag or more 
(|Herbielll977l ). Currently about 20 objects have been clas - 
sified as FU Ori- type (for a list see [Abraham et al.ll20o4 ). 
Records of the outbursts are not available for all of these 
objects: some were classified as FUors since they share many 
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spectral properties with the FUor prototypes. According 
to th e most widely accepted model IjKenvon fc HartmannI 
1 199ll) , the observed spectral energy distribution (SED) is 
reproduced by a combination of an accretion disc and an 
envelope. The energy of the outburst originates from the 
dramatic inc r ease o f the accretion rate. On the other hand, 
iHerbig et al] l|2003l ) favour a model in which the outburst 
occurs in an unstable young star rotating near breakup ve- 
locity. 

Parsamian 21 is a system consisting of a central ob- 
ject, HBC687 (^2000 = 19'' 29"" 0^87, (52ooo = 9° 38' 42'.' 69), 
and an ex t ended nebula first listed in the catalogue of 
IParsamianI (|l965l ). Though no optical outburst was ever 
observed, Parsamian 21 was classified as a FUor on the 
basis of its optical spectros copic and far-i nfrared proper- 
ties (Staude fc Neckell [l993 ). According to iHenning et al.l 
jlQQi ) Parsamian 21 is situated in a molecular cloud named 
"Cloud A", very close to the Galactic plane. Indeed, the dis- 
tance of 400 pc and the radial velocity of visr = -1-27 ± 
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15 km s ^ jStaude fc Neckell Il992l ') agree very well with 
the respective values for Cloud A (visr = +27 rfc 4 km 
s"\ rf = 500 ± 1 00 PC. iDame fc Thaddeud Il985h . Re- 
cently iQuanz et al.l (|2007l ) questioned the FUor nature of 
Parsamian 21 referring to mid-infrared spectral properties, 
which better resemble those of post-AGB stars (this issue 
will be discussed in Sect.|4]). 

The elongated nebulosity around the central source ex- 
tends ~ 1 arcmin to the north, while it is less developed 
to the south. It is visible in the optical and the near- 
infrared, though the source becomes unresolved at 10.8 
and 18.2 um (Poloms ki et alJ 120051 ). Polarimetric observa- 
tions by (Praper et al.l (|l985h revealed a centrosymmetric 



polarization pattern and an elongated band of low po- 
larization perpend i cular to the axis of the bright nebula. 
iBastien fc Menaxdl l|l990h interpreted the polarization map 
of Parsamian 21 in terms of multiple scattering in flattened, 
optically thick structures and derived an inclination angle 
of 80-85° and a size of 30 x 8 arcsec for this disc-like struc- 
ture. The existence of a disc is also supported by the discov- 
ery of a short bipolar o utflow oriented along the polar axis 
of the nebula jStaude fc Ncckcl 199^). Emissio n at su bmm 
wavelengths w a s det ected by iHenning et al.l (|l998l ) and 
IPolomski et~al] (j2005l l. The surroundings of Parsamian 21 
have been searched for companions several times, but close- 
by sta rs seen in J, H and K-band images proved to be field 
stars l|Li et al.l Il994l ). and no clos e-by sources have be en 
found at longer wavelengths so far (jPolomski et al.ll2005l '). 

Recently several studies were published on the cir- 
cumstellar environment of FUors using high-r esolution in- 
frare d and interferomet ric techniques (e.g. Green et al.l 



20061: iMalbet et al.l 1 19981 , l2005l : iMillan-Gabet et all l2006a 
Quanz et al.ll2007l ). Such observations could prove or dis 



prove the basic assumptions of FUor models. In this paper 
we present the highest resolution imaging and polarimetric 
data yet on Parsamian 21 from the VLT/NACO instrument, 
allowing the inspection of the circumstellar material and disc 
at spatial scales of ~ 0.07 arcsec. The edge-on geometry of 
Parsamian 21 represent an ideal configuration to separate 
the components of the circumstellar environment (disc, en- 
velope) and understand their role in the FUor phenomenon. 



2 OBSERVATIONS AND DATA REDUCTION 

The log of observations presented in this paper is summa- 
rized in Table [1] In the following subsections we describe 
these measurements in detail. 



2.1 VLT/NACO observations 

We have acquired ground-based near-infrared imaging and 
polarimetric observations in visitor mode using the NACO 
instrument mounted on the UT4 of ESO's Very Large 
Telescope VLT at Cerro Paranal, Chile on 2004 June 
18. NACO consists of the NAGS ada ptive optics system 
and the CON ICA ncar-infrarcd camera l|Lenzen et al.ll 19981 : 
iHartung et al. 2000; Roussct et al. 20031 ). The weather con- 
ditions were excellent, the typical optical seeing was ~ 0.72 
arcsec with two short peaks of 1.3 arcsec during the night. 
For imaging we used the visual dichroic and a camera with 
13 mas/pixel scale to obtain H (Ac = 1.66 /im) and Kg 



(Ac = 2.18 ^m) images and a camera with 27 mas/pixel 
scale to obtain L' (Ac = 3.8 /im) images. In each filter a 
4-point dithering was applied with small-amplitude random 
jitter at each location. The adaptive optics configuration was 
fine-tuned for each filter to ensure the best possible correc- 
tion. The H-band polarimetric observations were obtained 
using the 27 mas/pixel scale camera. 

VLT/NACO imaging. The data reduction was carried 
out using self-developed IDL routines. For the H and Kg 
band filters lampfiats were taken, while for the L' filter, 
skyfiats were acquired. For each filter, these images were 
combined to a final flatfield, and bad pixel maps were pro- 
duced. Then the raw images were flatfield corrected, and bad 
pixels were removed. Sky frames were calculated by taking 
the median of all images taken with the same filter, then im- 
ages were sky subtracted. Individual frames were shifted to 
the same position and their median was taken to obtain the 
final mosaic. Shifts were calculated by computing the cross- 
correlation of the images, which gives a more precise result 
in case of extended sources, than a simple Gaussian-fitting to 
the peak. Dithering resulted in a final mosaic of 21.4 x 21.4 
arcsec in case of H and Kg band, and 43.8 x 43.8 arcsec in 
L' band. The photometric standard was S889-E for the H 
and Kg filters (H = 11.662 ± 0.004 mag. Kg = 11.585 ± 0.005 
mag, IPersson et alii 1994)- and HD 205772 f or the L' fiher 
(L' = 7.636 ± 0.027 mag. lBouchet et al.lll99ll ). 

VLT/NACO polarimetry. Polarimetric observations 
were obtained using the differential polarimetric imaging 
technique (DPI, see e.g. [Draper et al.l Il985l : iKuhn et al.l 
I2OOII : lApai et al.ll200^ ). The basic idea of this method is to 
take the difference of two orthogonally polarized, simultane- 
ously acquired images of the same object in order to remove 
all non-polarized light. As the non-polarized light mainly 
comes from the central star, after subtraction only the po- 
larized light, such as the scattered light from the circum- 
stellar material remains. This standard way of dual-beam 
imaging polarimetry combined with adaptive optics at the 
VLT is a powerful way of doing high-contrast mapping. We 
obtained polarimetric images with NACO through the H fil- 
ter, using a WoUaston prism with a 2 arcsec WoUaston mask 
to exclude overlapping beams of orthogonal polarization. 
Parsamian 21 was observed at four different rotator angles of 
0°, 45°, 90° and 135°, providing a redundant sampling of the 
polarization vectors. At each angle a 3-point dithering was 
applied. The polarim etric calibrator was R CrA DC No. 71 
IWhittet et al.lll992l ). The polarimetric data reduction was 
done in IDL using previ ously developed so ftware tools which 
we presented in detail in lApai et al.l l|2004l ). According to the 
NACO User's Manual, instrumental polarization is generally 
about 2%. Based on our polarimetric calibration measure- 
ments, we found a typical instrumental polarization value of 
3% (position angle: 91° east of north). Considering the ex- 
pected high polarization for the Parsamian 21 nebula, we can 
conclude that instrumental polarization can be neglected in 
the present case, except in limited areas of low polarization. 

2.2 GFP Hc« imagery 

Parsamian 21 was observed with the Goddard Fabry-Perot 
(GFP) interferometer at the Apache Point Observatory 3.5m 
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Table 1. Log of the observations of Parsamian 21. Top part: data from this paper; Bottom part: archival data. 



Instrument 


Filter (central A) 


Mode 


Date 


Exp. Time 


Field of View 


Pixel scale 


FWHM 








YY/MM/DD 




[arcsec X arcsec] 


[arcsec] 


[arcsec] 


VLT/NACO 


H (1.66 /im) 


Imaging 


04/06/18 


8x30 s 


13.6x13.6 


0.013 


0.07 


VLT/NACO 


Kg (2.18 Atm) 


Imaging 


04/06/18 


8x20 s 


13.6x13.6 


0.013 


0.07 


VLT /NACO 


L' (3.80 /im) 


TmaerinET 


04/06/18 


48x0.2 s 


27.8x27.8 


0.027 


0.11 


VLT/NACO 


H (1.66 /im) 


Polarimetry 


04/06/18 


72x10 s, 24x80 


s 27.8x3.1 


0.027 


0.12 


GFP 


6563 and 6590A 


Imaging 


07/05/15 


900 s 


130 circular 


0.38 


1.1 


HST/NICMOS 


2 lira 


Polarimetry 


97/11/12 


18x31 s 


19.5x19.3 


0.076 


0.16 


HST/WFPC2 


F814W (0.80 Atm) 


Imaging 


01/07/30 


2x500 s 


79.5x79.5 


0.1 


0.17 


Spitzer/IRAC 


3.6, 4.5, 5.8 and 8.0Atm 


Imaging 


04/04/21 


848 s 


318x318 


1.2 


1.9-2.7 


Spitzer/MIPS 


24 and 70 /im 


Imaging 


04/04/13 


73 s, 462 s 


444x480, 174x192 


2.5, 4.0 


6, 18 


Spitzer/IRS 


5 — 40 lira 


Spectroscopy 


04/04/18 


1121 s 








ISO/ISOPHOT 


65 and 100 /xm 


Imaging 


96/09/28 


3542 s 


360x420 


15 


44, 47 



telescope on 2007 May 15. The observations were made in 
direct imaging mode. The on-band image has a central wave- 
length of 6563 A (He) and a width of 120 km s"\ The off- 
band image has a central wavelength of 6590 A and a width 
of 120 km s~^. The pixel scale is 0.38 arcsec per pixel. The 
observations were made through patchy cirrus and are not 
photometric, with seeing near 1 a rcsec. The inst r umen t and 
data reduction are described in IWassell et al. I l|2006l ) and 
references therein. For the Hq, on-band image, conspicuous 
night sky rings contaminate the region near Par 21. In order 
to remove these rings, we created a data cube consisting of 
subimages centred on the GFP optical axis, rotated by dif- 
ferent angles. After excluding the area around Parsamian 21, 
we azimuthally medianed the subimages, and subtracted the 
resulting image from the original measurement. For the off- 
band image, the night sky rings had displaced beyond the 
location of the Par 21 nebulosity, thus no ring subtraction 
was necessary. The on- and off-band images were shifted and 
scaled to match each other by measuring the positions and 
fluxes of seven stars in the vicinity of Parsamian 21. The 
continuum-subtracted Ha image can be seen in Fig. (3] 



2.3 HST archival data 

HST/NICMOS polarimetric observations were obtained 
on 1997 November 12 using the NIC2 camera in MULTI- 
ACCUM mode, and the POLOL, POL120L and POL240L 
polarizing filters. These filters have a bandpass between 1.9 
and 2.1 /im. Raw data files were calibrated at the STScI 
with the calmca v. 4. 1.1 pipeline. Since these data have not 
been published yet, we downloaded the pipeline-processed 
files and we did further processing using the IDL-based Po- 
larizer Data Analysis Software, which is available through 
the S TScI website and is described by iMazzuca fc HinesI 
l|l999 l) . This software package combines the images obtained 
th rough the three po larizers using an algorithm described 
in iHines et al.l (|2000l ). We used coefficients appropr i ate fo r 
the pre-NCS measurements from iHines fc Schneiderl (|2006D . 
A 6-point dither pattern was applied, covering in all about 
24 X 34 arcsec. The central star and the inner parts of the 



Parsamian 21 nebula can only be seen in two of the six im- 
ages, which we shifted and co-added. 

HST/WFPC2. Parsamian 21 was observed with the 
WFPC2 on 2001 July 30 through the F814W filter. 
Parsamian 21 is in the middle of the image on the WF3 chip 
(Fig. [T|, which has a pixel scale of 0.1 arcsec. Since these 
data have not been pubhshed yet, we used the High-Level 
Science Product with the association name 'U6FC2801B', 
available at the HST archive. 

2.4 Spitzer archival data 

Observations of Parsamian 21 were obtained on 2004 April 
13 (MIPS), April 18 (IRS) and April 21 (IRAC). The MIPS 
and IRAC, as well as part of the IRS data are still unpub- 
lished. The IRAC images cover an area of 5.3 x 5.3 arcmin 
centred on Parsamian 21. The 'high-dynamic-range' mode 
was used to obtain 15 frames in each position, 5 with 0.4 s 
exposure time and 10 with 10.4 s. The frames were processed 
with the SSC IRAC Pipeline vl4.0 to Basic Calibrated Data 
(BCD) level. The BCD frames were then processed using the 
IRAC artifact mitigation software, and finally the artifact- 
corrected frames were mosaicked with MOPEX. Photometry 
for Parsamian 21 was done in IDL. We used the short expo- 
sure images, because Parsamian 21 did not saturate these 
frames. We used an aperture of 6 arcsec and a sky annu- 
lus between 6 and 12 arcsec. The aperture corrections were 
1.061, 1.064, 1.067 and 1.089 for the four channels, respec- 
tiveljij. Colour correction was applied by convolving the ob- 
served SED with the IRAC filter profiles in an iterative way. 
The resulting colour-corrected fluxes can be seen in Table [S] 
Photometry for field stars was obtained in the long 
exposure images, because most field stars are fainter than 
Parsamian 21 and did not saturate even the long exposure 
frames. For this purpose we used StarFinder, which is an 
IDL-GUI based pro gram for crowded stellar fields analysis 
ijDiolaiti et al.|[2000l ). PSF-photometry was calculated for all 

^ IRAC Data Handbook, version 3.0, available at 
|http://ssc.spitzer. caltech.edu/irac/dh/ir acdatahandbookS.O.pdf I 
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stars having S /N > 3 and aperture corrections correspond- 
ing to "infinite" aperture radius were applied (0.944, 0.937, 
0.772 and 0.737 for the four channels, respectively^). Then, 
we cross-identified the sources in the four bands and found 
that 100 field stars are present in all four IRAC images. We 
plotted these sources as well as Parsamian 21 on a [3.6] — [4.5] 
vs. [5.8] — [8.0] colour-colour diagram in Fig. 1101 This plot will 
be discussed in Sect. 14.21 

The MIPS 24 and 70 /im images were taken in 'compact 
source super resolution' mode. At 24 /im 14 frames (each 
with 2.62 s exposure time), at 70 /im 44 frames (each with 
10.49 s exposure time) were combined into a final mosaic 
using MOPEX. At 24 /im Parsamian 21 saturated the de- 
tector, thus we used a model PSF to determine which pixels 
are still in the linear regime. Then we fitted the PSF only 
to these pixels. The 70 /^m image is not saturated, therefore 
we simply calculated aperture photometry in IDL using an 
aperture radius of 16 arcsec, sky annulus between 39 and 65 
arcsec, and an aperture correction of 1.7410. We found that 
the source was point-like at both 24 and 70 /^m. Colour cor- 
rection was applied by convolving the observed SED with 
the MIPS filter profiles in an iterative way. The resulting 
colour-corrected fluxes can be seen in Tabled 

The IRS spectroscopy was carried out in 'staring mode' 
using the Short Low, the Short High and the Long High 
modules. Measurements were taken at two different nod po- 
sitions and at each nod position three exposures were taken. 
We started from BCD level and extracted spectra from the 
2D dispersed images using the Spitzer IRS Custom Extrac- 
tion software (SPICE). In the case of the Short Low channel 
we extracted a spectrum from a wavelength dependent, ta- 
pered aperture around the star, and also from a sky position. 
In the case of the high-resolution channels, we extracted 
spectra from the full slit. For the Short High channel, the 
target was off-slit and most of the stellar PSF fell outside 
of the slit. We corrected this spectrum for flux loss using 
the measured IRS beam profiles (for more details, see the 
Appendix) . The resulting complete 5 — 35 /im spectrum can 
be seen in Figs. [S] and |S] 



Table 2. Photometry for Parsamian 21. All fluxes are colour- 
corrected. 



Instrument 


Wavelength [/tm] 


Flux [mjy] 


Spitzer/IRAC 


3.6 


143 ± 8 


Spitzer/IRAC 


4.5 


187 ± 10 


Spitzer/IRAC 


5.8 


303 ± 16 


Spitzer/IRAC 


8.0 


829 ± 42 


Spitzer/MIPS 


24 


5 530 ± 310 


Spitzer/MIPS 


70 


13 300 ± 1 300 


ISO/ISOPHOT 


65 


12 600 ±1300 


ISO/ISOPHOT 


100 


14 200 ± 1 900 




is^as^r 2" r or 



S8?=l5Sr 



Figure 1. HST/WFPC2 image of Parsamian 21 taken through 
the F814W filter. The intensity scale is square root and brightness 
increases from light to dark. The white square marks the area 
shown in the upper left panel of Fig. [2] 



2.5 ISO archival data 

ISO measurements of Parsamian 21 at 65 and 100 /tm 
were obtained in 'oversampled map' (PHT32) mode on 
1996 September 28. These observations were reduced us- 
in g a dedicated so f tware package (P32TOOLS) developed 
bv lTuffs fc Gabriell (|2003h . This tool provides adequate cor- 
rection for transients in PHT32 measurements. Absolute cal- 
ibration was done by comparing the source flux with the on- 
board fine calibration source. In order to extract the flux of 
our target, a point spread function centred on Parsamian 21 
was fitted to the brightness distribution on each map. At 
65 ^lm the PSF could be weU fitted by the ISOPHOT the- 
oretical footprint function. At 100 /im, however, the source 
turned out to be extended. In this case the brightness dis- 
tribution was modelled as the convolution of the standard 
theoretical footprint function with a two-dimensional ellip- 
tical Gaussian of 40 arcsec x 20 arcsec, and a position angle 



^ MIPS Data Handbook, version 3.2, available at 
|http://ssc. spitzer. caltech.edu/mips/dh/mipsdatahandbook3. 2.pdf I 



of 114°. The obtained fluxes were colour corrected by con- 
volving the observed SED with the ISOPHOT filter profiles 
in an iterative way. The results are displayed in Table [2] 



3 RESULTS 

3.1 Broad-band imaging 

Figure [U displays the HST/WFPC2 image taken at 0.8 /im. 
The image reveals the structure of the Parsamian 21 neb- 
ula with an unprecedented spatial resolution and detail. We 
note that strong emission lines associated with the bipolar 
outflow (see Sect. 13. 2p are not included in the bandpass, so 
the HST image provides a clean picture of the reflection neb- 
ula. The overall dimensions of the nebula are approximately 
60 X 20 arcsec. The northern part has a characteristic el- 
liptic, loop-like shape, with the inside of the loop having a 
relatively low surface brightness. The outer edge of the neb- 
ula is rather well-defined, while the inner edge is more fuzzy. 
Apart from the loop itself, there is also a wide, faint SE-NW 
oriented stripe about 30 arcsec from the star. The nebula is 
clearly bipolar, although very asymmetric: the southern part 
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Figure 3. Continuum-subtracted Hq image of Parsamian 21. The 
peak at position (0,0) is the residuum of the central source. Three 
Herbig-Haro knots could be identified; these were named Knot 1, 
2 and 3 bv iStaude fc Neckcl 1,1992). 



Figure 2. HST/WFPC2 (with filter F814W) and VLT/NACO 
(with filters H, Kg and L') images of Parsamian 21. Circles at 
the left bottom corners indicate the FWHM at the centre of the 
corresponding image. In the Kg-band image, arrows mark the 
positions of the two closest stars to Parsamian 21 (see Sect. l4?2t . 

is much less developed than the northern one. We note that 
while the southeastern arc is brighter close to the star (see 
the central part of the HST image in Fig. [2| , the southwest- 
ern arc is more conspicuous farther from the star (Fig. [1]). 
It is remarkable that field stars are visible both to the north 
and to the south of Parsamian 21. Supposing that these stars 
are in the background, this implies that there is no signifi- 
cant difference in the extinction between the northern and 
southern part. 

Figure [2] shows the central parts of the HST/WFPC2 
and the VLT/NACO images (note that the innermost 0.57 
arcsec of the HST image is saturated). The bipolar nature of 
the nebula is most conspicuous at the shortest wavelength 
(0.8/.im), where all 4 arcs (roughly to the northeast, north- 
west, southwest and southeast) are visible. In the Ks and H 
bands the nebula is more triangle-shaped. In all three NACO 
bands the central source is extended, with deconvolved sizes 
of 0.20, 0.12 and 0.08 arcsec in the H, Ks and L' bands, 
respectively. 

3.2 Narrow-band imaging 

Figure |3] displays the continuum-subtracted image of 
Parsamian 21. Conspicuous Herbig-Haro knots are seen to 
the north and the south of the central star. The north- 
ern knot is compact, resembling the knot described in 
IStaude fc Neckell (| 19921 ). The southern knot, however, has 
a distinctly bowed shape. The southern extent of this knot 
has a S/N of 10 over background per pixel, while the fiank- 
ing structures have S/N ~ 7. The northern knot has S/N 
in excess of 15 per pixel. To measure the angular distance 
of the knots from the central star in the H CK imagery, we 
binned the data in a 3 pixel (1.1 arcsec) wide swath ori- 
ented north-south. Apart from the residual of the central 



1200 - 




-10 -5 5 

Angular distance [arcsecl 



Figure 4. South-north oriented cut across Parsamian 21 in the 
continuum-subtracted Hq image. marks the position of the cen- 
tral star. The peaks marked by solid vertical lines at —10.5, —3 
and 5 arcsec correspond to Knot 2, 3 and 1, respectively. Dashed 
vertical lines mark the p ositions of the same knots as measured by 
IStaude &: Necke] lll992l ) . Arrows demonstrate the direction where 
the knots propagate. 

source, we could identify three knots (see Fig. [S}: Knots 1, 

2 and 3 are located 5 arcsec north, 10.5 arcsec south and 

3 arcsec sout h, respectively. The knot s correspond to those 
identified by IStaude fc Neckell l| 19921 ). but all three knots 
have moved outwards since then (Fig. |3| . Assuming a dis- 
tance of 400 pc, and supposing that the projected velocities 
are close to the velocities of the knots along the polar axis, 
one can derive 120 km s~^ for Knots 1 and 3, and 530 km 
s~^ for Knot 2. This is in good agreem ent with the velocities 
calculated by IStaude fc Neckell (ll992') using radial velocity 
measurements. Such outfiow velocities imply kinematic ages 
of approximately 80 yr for Knot 1 and 40 yr for Knots 2 and 
3. 

3.3 Spectral energy distribution 

Table [2] contains our Spitzer and ISO photometry for 
Parsamian 21. These data, complemented with previously 
published measurements between 1983 and 1999 are plot- 
ted in Fig. O showing the complete UV-to-mm spectral en- 
ergy distribution (SED) of the object. Comparing data taken 
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Figure 5. Complete UV-to-mm SED of Parsamian 21. Source of 
data: 2MASS All-Sky C atalog of Point S o urces , MSX6C Infrared 
Point Source Catalog, [A braham et al.l ||200' 
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||2005| '). iNeckel fc Stand J (1984). Hcnni ng et"al1 {l998,) and this 
work. The model SED is discussed in Sect. 14.41 
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Figure 6. 5 — 18 /.tm part of the Spitzer/IRS spectrum of 
Parsamian 21. The inset displays in details the continuum- 
subtracted spectral region around the 15.2 /im CO2 ice feature 
of Parsamian 21 with solid line and as a comparison HH46IRS 
with dashed line (for discussion see Sect. l43t . 



between 1983 and 1999, lAbraham et all (|2004l ') found that 
there is no long-term flux variation in the 1 — 100 j-im wave- 
length range. The new Spitzer data from 2004 reveal that 
the mid-infrared brightness of Parsa mian 21 stayed constant 
since then at a 15% level. Although IParsamian et al.1 (|l996l ) 
reported 2 — 3 mag brightness variations in the B band be- 
tween 1966 and 1990 (in Fig. [5] we plotted optical measure- 
ments from 1980 as representative values) , no significant flux 
changes can be seen at longer wavelengths within the mea- 
surement uncertainties. 

At w avelengths shorter th an 3 fim, the SED is highly 
reddened l|Polomski et al.ll2005l ). probably due to the com- 
bined effect of interstellar extinction and self-shadowing by 
circumstellar material. The fact that the central source is 
extended at these wavelengths (see Sect. I3.l"|l indicates that 
the photometry is contaminated by scattered light from 
the inner part of the circumstellar environment. Between 
3 and 25 ^m, the SED is rising towards longer wavelengths 
as uFi, oc A. The 5 — 18 ^m Spitzer spectrum (Fig. [6} dis- 
plays strong PAH emission features at 6.3, 8.2 and 11.3 /xm 



and a morphous silicate emissi on around 10 ^m (("Quanz et al.l 
|2007| . see also the spectrum of lPolomski et al.ll2005h . Weaker 
absorption features can also be seen between 7 and 8 fim, 
which cannot be safely identified with any ice or PAH 
feature. The 10 — 19 // m channel of t he Sp itzer/IRS spec- 
trum was not shown bv lQuanz et al.l l|2007l ) due to data re- 
duction difficulties. We made an attempt to evaluate also 
this channel (for details, see Appendix). As can be seen 
in the inset of Fig. [S] the spectrum shows a strong CO2 
ice absorption feature at 15.2 fim, which has not been re- 
ported in the literature so far. Between 25 and 100 fim the 
SED is flat (vF^ cx const.), while the submillimetre shape 
is uFiy cx A~^'^ (corresponding to /3 = 0.5, assuming optically 
thin emission and a dust opacity law of k,/ oc i^*^). Using a dis- 
tance of 400 pc and an in terstellar reddening of Av ~ 2 mag 
iHiUenbrand et all [l993 ). the bolometric luminosity com- 
puted as the integral of the SED from 0.44 to 1300 fim is 
10 L0 (assuming isotropic radiation fleld). 



3.4 Polar imetry 

Figure [7] shows images of Parsamian 21 calculated from the 
H-band NACO polarization data, as well as from the 2 /im 
NICMOS measurements. The left column shows the total 
intensity (/) obtained as the sum of the intensities of the two 
orthogonal polarization states. The middle column displays 
the polarized intensity (PI) calculated as the square root 
of the quadratic sum of the Q and U Stokes components. 
The right column shows the degree of polarization (% poL), 
i.e. the ratio of the polarized to the total intensity. 

The NACO total intensity maps (Fig. [T] left) are very 
similar to the H-band map in Fig. (2] Thus, the WoUaston- 
prism data can reproduce well the direct images taken with- 
out the prism. While the total intensity is relatively smooth, 
the polarized intensity (Fig. [71 middle) reveals many fine de- 
tails not visible in the total intensity map. The most strik- 
ing feature is a dark horizontal lane across the star, where 
the polarized intensity is very low. Above this dark lane, 
high intensity regions can be seen, primarily emphasizing 
the walls of the upper lobe. Below the dark lane, the south- 
eastern arc is very pronounced. The degree of polarization 
(Fig. [3 right) again shows features different from the previ- 
ous two maps. Here, the central star almost disappears, as 
the stellar light is not polarized. The horizontal dark lane 
of low polarization across the star is even more pronounced: 
there the degree of polarization is around ~ 5% in the H- 
band NACO images. We note that the polarization pattern 
in this area is probably affected by instrumental polarization 
(Sec. 12. l) . though this effect should not alter our discussion 
in Sec. 14.3.21 The polarization of the central source itself, 
measured in a 0.08 arcsec radius aperture, is 9% and its po- 
sition angle is 70° east of north. The long exposure NACO 
images (Fig. [3 middle row), showing the outer regions with 
higher signal-to-noise ratio, are consistent with the short ex- 
posure maps (Fig. [T] upper row). It is remarkable that the 
NICMOS images, which represent independent observations 
(different resolution, different instrument, different polariza- 
tion technique), show strikingly similar features, although 
the degree of polarization is slightly higher. 

Figure |S] shows the polarization pseudovectors overlaid 
on total intensity contours. The NACO map, along with the 
degree of polarization map in Fig. [7] right shows that the 
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Figure 7. Left: total intensity (/); middle: polarized intensity (PI); right: degree of polarization (% pol.), with scalebar. Upper row: 
short exposure H-band NACO images; middle row: long exposure H-band NACO images (the saturated central parts are masked out); 
lower row: 2 /xm NICMOS images of the same area. The intensity scale is logarithmic and brightness increases from dark to light. 



horizontal band of low polarization across the star is very 
well-confined, narrow (about 0.15 — 0.20 arcsec wide at the 
5% polarization level) and can be discerned from 0.12 to 
1.2 arcsec to the star (and likely beyond, as indicated by 
the long exposure ma ps, see also the po larization maps of 
[Draper et al]ll985l and lHaiiar et al.lll997l '). The position an- 
gle of this band i s 78° ± 4° east of n orth, consistent with 
that measured bv lDraper et all l|l985h (75° ±4°). This im- 
plies that the band in our high spatial resolutio n images are a 
direct inwa rd continuation of t he band seen bv lDraper et al.l 
l|l985l ) and lHaiiar eFal] (119971 1 in their lower resolution im- 
ages. The degree of polarization in this band is a few percent, 
and it is oriented parallel with the band itself. 

Most of the northern part of the nebula shows a 
centrosymmetric pattern, characteristic of reflection neb- 



ulae with single sc atteri ng. This is consist ent with what 
iDraper" et all l|l985l ) and iHaiiar et al.1 (|l997l ') measured for 
Parsamian 21 itself. It is also si milar to what can b e seen in 
other bipolar nebulae (see e.g. 'Mcakin et al.l [20051 for neb- 
ulae of young stars or Scarrott ct al. 199^ for nebulae of 
evolved stars). The highest polarization in the northern part 
can be measured in the northeastern wall (25 — 30%) and 
also in a spot about 0.9 arcsec west and 0.6 arcsec north 
to the star (20 — 25%). The southern part of the nebula 
also follows the regular centrosymmetric pattern seen in the 
northern part. Here the highest polarization can be found 
in the southeastern arc (60 — 70%). 

The NICMOS vector map at 2 /im (Fig. [S] bottom) re- 
peats all the main features seen by NACO in the H-band 
(Fig.[8]top). We note nevertheless the presence of two small. 
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Figure 8. Polarization pseudovectors overlaid on intensity con- 
tours. Top: VLT/NACO, H-band, with 5x5 pixels binning, bot- 
tom: HST/NICMOS, 2/im, with 2x2 pixels binning. 



localized depolarization areas on either side of the central 
source at positions ((/.'g, (/.'O) and (-(/.'S, -O'.'l). 



4 DISCUSSION 

4.1 Parsamian 21: an FU Orionis-type star 

Although no optical outburst was ever observed, 
Parsamian 21 shares many properties characteristic of 
FUors. Its spectral typ e in the liter ature ranges from A5 
to F8 superg iant (e.g. IStaude fc Ne ckcl 1992), similar to 
the typical FUor spectral types (F-G supergiant). In ad- 
dition, it shows strong infrared excess and drives a bipolar 
outflow, similarly to many FUors. However, in a recent 
paper bv lQuanz et al.l (|2007l ) the pre-main sequence nature 
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Figure 9. 2MASS colour-colour diagram of 188 sources in an area 
of 5.3 X 5.3 arcmin centred on Parsamian 2 1. Th e main-sequence 
and giant branch are marked by sol id lines llKoo rnnccf 1983), the 
reddening path with dashed li nes llC ardelli etal.i ll989[') and the 
TTauri locus with dotted line l lMever et al.lll997t) . 



and the FUor status of Parsamian 21 were questioned, 
mainly because its PAH emission bands are untypical for 
young stars. They suggest that Parsamian 21 is either an 
intermediate mass FUor object, or an evolved star sharing 
typical properties with post-AGB stars. 

An important parameter which may discriminate be- 
tween a pre-main sequence object and a post-AGB star is 
the luminosity. The typical lum inosity of p ost-AGB stars 
is in the order of 10^ - 10* Lq (|Kwoklll993l ). while FUors 
are typically fainter than a few 100 L©. The luminosity of 
Parsamian 21, 10 L0 (see Sect. 13. 3p . is at least a factor 
of 100 lower than that of post-AGB stars. Even adopt- 
ing the largest distance estimate men tioned in the liter- 
ature (1800 pc, IStaude fc Neckell liggj ). the luminosity of 
Parsamian 21 is still too low. The observed proper motion of 
the Herbig-Haro knots, however, prefers the lower distance 
(400 pc), otherwise the outflow velocities (« 2400 km s~^ at 
1800 pc) would be unusually high for a young stellar object 
(typically up to 600 km s"\ see e.g. iMundt et"al] Il987l : 
iHessman et al.l Il99ll) . Moreover, the existence of Herbig- 
Haro outflows are usuall y tracers of low-mass star formation 
iReipurth fc Ballvll200ll ). 

In the following discussion, we consider Parsamian 21 
to be a FUor and we discuss its properties in the context of 
young eruptive stars. Nevertheless we note that most of our 
results on the morphology and circumstellar structure are 
valid regardless of the nature of the central object. 



4.2 Parsamian 21: an isolated young star? 

Parsamian 21 is situated close to the Galactic plane (I — 
45.8°, 6 — —3.8°), in a molecular cloud calle d Clou d A. This 
cloud was identified by iDame fc ThaddeusI (|l985D in their 
CO survey of molecular clouds in the northern Milky Way. 
The cloud occupies an area of 8 square degrees in the sky. 
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Figure 10. IRAC colour-colour diagram of 100 sources located 
in the same area as in Fig. |9] The gray square in the upper 
right corner mark s the approximate domain of Class II sources 
llAllen et al.ll2004h . 




Figure 11. Sketch of the morphology of circumstellar material 
around Parsamian 21, overlaid on the HST/WFPC2 image. The 
central star is surrounded by an edge-on disc. Perpendicular to the 
disc, the star drives a bipolar outflow that excavates an outflow 
cavity in the dense circumstellar material. Light from the central 
star illuminates the walls of the cavity. 



Table 3. VLT/NACO photometry for the two closest stars. Un- 
certainties are about 0.1 mag. 
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and the only known young star associated w ith it is the 
TTauri star AS 353 (|Dam6 fc Thaddeusll 19851 ). 

In order to check whether FUors are usually associated 
with star forming regions, we searched the literature and 
found that mos t FUors are located in areas of active star 
formation (e.g. iHenning et al.|[l998l ). To find out whether 
there is star formation in the vicinity of Parsamian 21, 
we searched for pre-main sequence stars. For this purpose 
we constructed a 2MASS J-H vs. H-Ks and an IRAC 
[3.6] — [4.5] vs. [5.8] — [8.0] colour-colour diagram for sources 
found in our 513 x 5.'3 IRAC field of view (Fig. |l[10l). Our 
selection criteria in the case of 2MASS was S/N > 10 and 
uncertainties < 0.1 mag in all J, H and K5 bands, while 
in the case of IRAC S/N >3 and detectability at all four 
bands were required. The 2MASS diagram revealed that 
most of the nearby objects are reddened main sequence or 
giant stars. On the IRAC diagram, however, there are three 
objects (apart from Parsamian 21 itself), which display in- 
frared excess at 8 /im (marked b y A, B and C in F ig. fUJ)) . 
According to the classification of I Allen et al.l l|2004l ). Class 
II sources exhibit colours of [3.6] -[4.5] > 0.0 and [5.8] -[8.0] 
> 0.4, thus one of these stars (A) might be a Class II source, 
while B and C are more likely Class Ill/main sequence 
sources. The nature of source A and its possible relation- 
ship to Parsamian 21 is yet to be investigated. Nevertheless, 
Parsamian 21 seems to be rather isolated compared to most 
FUors and certainly not associated with any rich cluster of 
young stellar objects. 

We also searched for possible close companions of 
Parsamian 21 in the WFPC2 and NACO direct images. In 
order to establish a detection limit for source detection, we 
measured the sky brightness on the NACO images (before 
sky-subtraction), and estimated a limiting magnitude for 



each filter. The resulting values are 22.8, 21.6, and 15.2 mag 
in H, Ks and L', respectively. In case of the HST/WFPC2 
image, the larger (80 x 80 arcsec) field of view made it pos- 
sible to estimate a limiting magnitude using star counts; 
the resulting value is 23.5 mag. Due to the bright reflec- 
tion nebula, the detection limit is somewhat lower close to 
the star. The two closest objects we found are the follow- 
ing: one star to the southeast, at a distance of 1.4 arcsec 
(560 AU at 400 pc), and another one to the northwest, at 
a distance of 3.3 arcsec (1320AU at 400pc). These sources 
are marked with arrows in Fig. [5] Neither of the stars are 
visible at 3.8 or 0.814 ^m, although we can give an upper 
limit for their L' brightness. Their positions and photom- 
etry are given in Table O As these sources are very red, 
they can equally be heavily reddened background stars, or 
stars with infrared excess (indicating that they might be as- 
sociated with Parsamian 21). Supposing that they are red- 
dened main sequence stars, one can estimate an extinction 
of Av ~ 10 — 15 mag. Further multifllter observations may 
help to clarify the nature of these objects and their possible 
relationship to Parsamian 21. 

4.3 The circumstellar environment of 
Parsamian 21 

The appearance and polarization properties of the nebula 
around Parsamian 21 can be understood in the following 
way: the star drives an approximately north-south oriented 
bipolar outflow, which had excavated a conical cavity in 
the dense circumstellar material (Fig. Ilip . The star illu- 
minates this cavity and the light is scattered towards us 
mainly from the walls of the cavity. The outflow direction 
is perpendicular to an almost edge-on dense circumstellar 
disc. This picture is supported by the following facts: (a) the 
centrosymmetric polarization pattern is characteristic of re- 
flection nebulae with single scattering; (b) the morphology 
and limb brightening suggest a hollow cavity (as opposed 
to an "outflow nebula", where the lobes are composed of 
dense material ejected by the central source); and (c) the 
low-polarization lane across the star strongly suggests the 
presence of an edge-on circumstellar disc, where multiple 



10 A. Kospdl et al. 




1 0.5 -0.5 -1 

RA offset ["] 



Figure 14. VLT/NACO polarization map of Parsamian 21 overlaid on H-band total intensity contours. The circle in the upper right 
corner displays the FWHM of the polarization measurement. Polarization vectors are displayed at full resolution, only showing the central 
low polarization band, where the polarization vectors are aligned. Such arrangement is expected when multiple scattering occurs in an 
edge-on disc. 
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Figure 12. Brightness profiles of Parsamian 21 at 0.8 ^m and in 
the H and Kg bands. Dotted lines mark Hubble's relation. 
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Figure 13. South-north cut at 0.6 arcsec east from the star. Solid 
line: total intensity; dashed line: polarized intensity. 



scattering occurs. In general the Parsamian 21 system shows 
similarities to the NGC2261 nebula associated with RMon. 
This object also consists o f a northern cometary nebu la and 
a southern jetlike feature (IWarren-Smith et al.lll987h . 

4-3.1 Envelope /Cavity 

We characterised the opening of the upper lobe by mark- 
ing the ridge along the northeastern and northwestern arcs 
which we interpret as the walls of the cavity. As viewed from 
the star northwards, the cavity starts as a cone with an open- 
ing angle of « 60°, giving the nebula in Fig. [5] a characteris- 
tic equilateral triangle-shape. Farther away from the star the 
cavity deviates from the conical shape, becomes narrower. 



The whole cavity occupies an area of 8 000 x 24 000 AU (at a 
distance of 400 pc). The sharp outer boundary of the nebula 
implies a significant density contrast between the cavity and 
the surrounding envelope. 

In Fig.[T2]we plotted radial brightness profiles at 0.8 /im 
and in the H and Ks bands, starting from the star north- 
wards. The slope of the intensity profiles in t he inner ~ 1.6 " 
(640 AU) follows closely Hubble's relation l|Hubblel [l92^ ) ■ 
i.e. the brightness is proportional to r~^. This trend can be 
clearly followed above the noise out to about 5 arcsec in our 
H and Ks images, giving an estimate of 2000 AU for the 
outer size of the envelope. The fact that the brightness pro- 
files follow Hubble's relation implies that the nebula is pro- 
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duced by isotropic single scattering, in accordance with the 
centrosymmetric polarization pattern and the high degree 
of polarization. Around ~ 1.6" the profiles becomes steeper, 
probably due to decreased density in the inside of the cav- 
ity. Similar steepening in the outer part of the nebula was 



Table 4. Model parameters. Parameters in italics are fixed, while 
the others were fitted. 



already mentioned bv lLi et al.l (|l994l ). The profiles are sim- 
ilar at all observed wavelengths and do not show significant 
dependence on the position angle. 

i.3.2 Disc 

As mentioned in Sect. [SJl both the NACO and NICMOS 
polarimetric images show a lane of low polarization oriented 
nearly east-west across the star (Fig. (Tj. The drop in the 
degree of polarization can also be clearly seen in Fig. 1131 
where a north -south cut at 0"6 e ast t o the star is plot- 
ted. Following iBastien fc MenardI (|l990l ). we interpret the 
low polarization by multiple scattering in an edge-on disc 
(possible other explanations for the origin of the low po - 
larization areas are discussed in e.g. iLucas fc RochdlTooi ). 
According to models of such circumstellar stru ctures (e.g. 
IWhitnev fc Hartmannlll993l : iFischer et al.lll996l ) the polar- 
ization vectors are oriented parallel to the disc plane. As can 
be seen in Fig. 1141 despite the low degree of polarization, the 
predicted alignment of the vectors can be clearly seen in the 
case of Parsamian 21 too. The NICMOS polarization map 
shows a similar effect. The dark lane in Fig. [7] can be fol- 
lowed inwards to as close as 48 AU from the star. This is an 
upper limit for the inner radius of the circumstellar disc. 

An interesting feature of the mentioned models is two 
depolarized areas on either side of the central star, which 
mark the outer end points of an edge-on disc. The depolar- 
ization is due to a transition from the linearly aligned to the 
centro-symmetic polarization pattern. These depolarized ar- 
eas can be seen in Fig.[8]bottom, on both sides at about 0.9" 
from the star. At the distance of Parsamian 21 this corre- 
sponds to 360 AU, and can be adopted as the outer radius 
of the dense part of the disc, where multiple scattering at 
near-infrared wavelengths is dominant. It is interesting that 
the NACO map does not show clear depolarized areas at the 
same position, but seems to resolve the transition in vector 
orientation, while the larger beam of NICMOS averaged the 
differently oriented vectors, resulting in depolarized spots. 

The polarized intensity and degree of polarization im- 
ages in Fig. [7] suggest a slight asymmetry in the dense disc: 
the eastern (left) side is straight, while the western (right) 
side shows a kink and also has a different position angle than 
that on the other side. The thickness of the disc can be mea- 
sured on the area where the polarization vectors are aligned 
(Fig. I14|l . This approximately corresponds to the area where 
the degree of polarization is below ~ 10%. The resulting 
thickness is approximately 0.1" (40 AU) to the east and is 
somewhat larger, 0.2" (80 AU), to the west. One should note 
that, since these values are close to the spatial resolution of 
the polarimetric images, these numbers should be consid- 
ered as upper limits for the thickness of the circumstellar 
disc around Parsamian 21. They are upper limits also be- 
cause if the inclination is not exactly 90°, the thickness of 
the disc can be even less. The thickness does not show signifi- 
cant increase with radial distance, suggesting the picture of a 
flat, rather than a flared disc, at least considering the dense, 
multiple-scattering part. The smooth brightness and polar- 
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Value 


Inner disc radius 


Ri 


3.5^0 


Outer disc radius 


R2 


360 A U 


Temperature at 1 AU 


Td.o 


285 K 


Power-law index for temperature 


<?d 


0.75 


Power-law index for surface density 


Pd 


1.6 


Disc mass 


Aid 


0.02 Mq 


Inclination 


i 


86° 


Inner envelope radius 


Rs 


5.4 AU 


Outer envelope radius 


Ra 


2000 A U 


Temperature at 5 AU 


Tc,o 


368 K 


Power-law index for temperature 


go 


0.4 


Power-law index for surface density 


Po 


0.4 


Envelope mass 




0.02 Mq 


Interstellar Extinction 


Av 


2 mag 



ization distribution between the disc and the surrounding 
envelope (Fig. I13p . however, implies that there is a contin- 
uous density transition between the two components. From 
the ratio of the horizontal to vertical sizes of the disc a lower 
limit of 84° for the inclination of the system (the angle be- 
tween the normal of the disc and the line of sight) can be 

deri ved. 

iFischer et al.l (|l996l ) computed a grid of polarization 
maps of young stellar objects with the aim of help- 
ing the interpretation of polarimetric imaging observa- 
tions. They consider five different models, four with mas- 
sive, self-gravitating discs and one with a massless Ke- 
plerian disc. Since the mass of circumstellar material of 
Parsamian 21 derived f rom submillimetre observations is rel- 
ative l y low (<0.3M(:), Henning et al. 1998!;_ Polomski et al.l 
I2OO5I : ISandell fc Weintraubll200ll : IHiUenbrand et al.lll992l ). 
the most appropriate model for our case is a Keple- 
rian disc. Ind eed the polarization pattern as computed by 
IFischer et al] ([l996i ') for an inclination of 87° (their Fig. 1) 
looks remarkably similar to our Fig. |S] (the differences might 
be explained by the narrower cavity and flatter disc of 
Pa rsamian 21). Thus, the geometry and structure assumed 
bv IFischer et al.l |l996) in their Keplerian model could be a 
good starting point for further radiative transfer modelling 
of Parsamian 21. 



4.4 Modelling the circumstellar environment 

Our observations provide some direct measurements of the 
geometry of the circumstellar structure (disc size and thick- 
ness, inclination, envelope size). In the following we discuss 
the consistency of this picture with the observed SED. Our 
approach is to construct a simple disc-|-envelope model, in 
which we flx those parameters whose values are known from 
our NACO observations or from other sources (outer disc 
radi us from this work; power-la w index for disc temperature 
from lShakura fc Sunvaevll 19731 : disc mass was set in order to 
ensure that the whole disc is optically thick; inclination from 
this work; outer envelope radius from this work, the power- 
law index for envelope temperature is a typical value for 
optically t hin envelopes co ntaining larger than interstellar 
grams, e.g. iHartmannlbOOd Eqn. 4.13; interstellar extinction 
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from iHillenbrand et al.lll992l ). Then we check whether the 
SED can be fitted by tuning the remaining paramet ers. We 
adopted an analytical disk model l) Adams et al.ll987[). which 
has been successfully used to model FUors ([Quanz et alj 
I2OO6I : [Abraham et alJliooil ) . Our model consists of two com- 
pone nts, an optically thick and Reometrically thin accretion 
disc (jShakura fc SunvaevI llOTSi ) and an optically thin en- 
velope (no cavity is assumed). No central star is included 
in the simulation, partly because in out bursting FUors the 
star's co ntribution is negligible com pared to that of the in- 
ner disc (iHartmann fc KenvonlHooil ) . and partly because of 
the edge-on geometry where the star is obscured by the disc. 
This assumption is supported by the fact that the shape of 
the SED at optical wavelengths is broader than a stellar pho- 
tosphere. The model also does not take into account internal 
extinction and light scattering, thus it cannot reproduce any 
of the near-IR imaging and polarimetric observations. 

The temperature and surface density distribution in the 
disc are described by power-laws: 



T(r) 



E(r) 



1 AU 



(liu) 



(1) 



(2) 



Similar power-laws were assumed for the envelope. The ob- 
served flux at a speciflc frequency is given by 

= ^ 27rr(l-e^^)B.(Td)dr+ (3) 



R4 



— / 27rr(l-e-^"-')B,(re)dr. 

^ JR3 
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The flrst term describes the emission of the accretion disc, 
the second term describes the radiation of the optically thin 
envelope. For the dust opacity we used a constant value 
of = 1 cm^g~^ at A > 1300 fim, = KiaooAim ( isoo^m ) 
between 1300 and 100 ^.m and again a constant value of 
= ft;(100 /im) at A < 100 ^.m. 
We fitted the SED via minimisation using a genetic 
optimization algorithm PIKAIA (jCharbonneau 1995'). This 
algorithm performs the maximization of a user defined func- 
tion, for which purpose we used the inverse x'^- Since there 
are many photometric measurements in the mid-infrared do- 
main, but just a few in the far-infrared, the mid-infrared 
region has a higher weight during the fit, compared to the 
far-infrared domain. Therefore, in order to ensure an equally 
good fit at all wavelengths, we divided the SED into four re- 
gions and weighted the of each domain with the inverse 
of number of photmetric points the region contained. Then 
the final was the sum of the x^ of all regions. The regions 
we used were: 0.3 - 3, 3 - 30, 30 - 300 and 300 - 3000 ^m. 
The parameters of the best-fit model, which gives a weighted 
X^ of 0.67, are fisted in Table H The fitted model SED as 
well as the disc and envelope components are overplotted in 
Fig. El 

The model SED is consistent with the observed fiuxes. 
This shows that the picture of a thin accretion disc and 
an envelope is consistent with both the measured SED and 
the geometry and disc/envelope parameters inferred from 
our polarimetric observations. Detailed modelling of the sil- 
icate, PAH, and ice spectral features, as well as the correct 
treatment of internal extinction and scattering would re- 



quire radiative transfer modelling, which will be the topic of 
a subsequent paper. 

4.5 The evolutionary status of Parsamian 21 

The geometry of FUor models discussed in the literature 
(e.g. Hartmar m &: Kenvo n 1996; Turner et al. 1997i) usually 
consist of a central star surrounded by an accretion disc 
and an infalling envelope with a wind-driven polar hole. 
These assumptions are suppor t ed by the fact that the y 
fit well the SED (ICreen et all I2OO6': 'Quanz et al.' '200 ' 
t he interferometric visibilities jMil lan-Gabct ct al. 2006_ 
Abr aham et al l2006l\ and the temporal evolution of the SED 
(|Abraham et al.l 20041 ') . In this paper we present the first di- 
rect imaging of these circumstellar structures in a FUor. Our 
polarimetric measurements of Parsamian 21 show the exis- 
tence of a circumstellar disc which extends from at least 48 
to 360 AU. The most striking feature of the disc is its flat- 
ness over the whole observed range. The short-wavelength 
part of the SED could be well reproduced using a radial 
temperature proflle of r""'^^ (Sect. [4^ . This proflle is ex- 
pected from both a geometrically thin accretion disc and a 
flat reprocessing disc. An envelope was also seen in the po- 
larization maps of Parsamian 21 and it was also a necessary 
component for the SED modelling. Envelopes are involved 
in many FUor models and in this paper we present a direct 
detection of this model component. Our images reveal that 
the envelope can be followed inwards as close to the star as 
the disc. FUor models often assume a polar cavity in the 
envelope, created by a strong outflow or disc wind. The di- 
rect images of Parsamian 21 clearly show the presence of 
such a cavity and we also detected a bipolar outflow in the 
Parsamian 21 system. 

In the recent years, as new interferometric and infrared 
spectroscopic observations were published for FUors, the 
group turned out to be more inhomogeneous in physical 
properties than earlier assumed, when mainly optical pho- 
t omet ry and spectroscopy had been available. lOuanz et al.l 
(200?) proposed that some differences might be understood 
as an evolutionary sequence. They suggest that FUors con- 
stitute the link between embedded Class I objects and the 
more evolved Class II objects. Members of the group exhibit- 
ing silicate absorption at 10 /xm are younger and more em- 
bedded (Category 1, e.g. V346Nor); while objects with pure 
silicate emission are more evolved (Category 2, e.g. FU Ori 
and Bran 76). There are objects showing a superposition of 
silicate absorption and emission, which are pro bably in an 
interm ediary evolutionary stage (e.g. RNO IB) . Green et al.l 
(|2006t l also sorted FUors, based on the ratio of the far- 
infrared excess and the luminosity of the central accre- 
tion disc, fd (Equ. 7 in their paper). A large relative ex- 
cess [fd > 5%) indicates an envelope of large covering frac- 
tion (V1057Cyg and V1515Cyg), while low relative excess 
means a tenuous or completely missing envelope (Bran 76 
and FU Ori). This is also an evolutionary sequence, as young, 
more embedded objects have large envelopes, while around 
more evolved stars, the envelope has already dispersed, fd 
can also be used to calculate the opening angle of the en- 
velope, thus a prediction of this scheme is that the opening 
angle is becoming wider during the evolution, probably due 
to strong outflows during the repeated FUor outbursts. 

The two classification schemes are not inconsistent and 
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one can merge them into the following evolutionary se- 
quence: (1) the youngest objects exhibit silicate absorp- 
tion and large far-infrared excess (V346Nor, probably also 
OOSer and L1551IRS 5 belong here); (2) intermediate- aged 
objects, where the silicate feature is already in emission but 
there is still a significant far-infrared excess (V1057 Cyg, 
V1515Cyg, probably also RNO IB and V1647 0ri); (3) the 
most evolved objects show pure silicate emission and low far- 
infrared excess (FUOri, Bran 76). We note, however, that 
this c lassification has some weak points. As iQuanz et al.l 
|2003) already mentioned, an edge-on geometry in a more 
evolved system may appear as a younger one. Moreover, 
during an outburst and the subsequent fading phase, cer- 
tain spectral features as well as the global shape of the SED 
may change. 

Parsamian 21 can be placed in this evolutionary scheme, 
though one should keep in mind that because of the nearly 
edge-on geometry, the classification of this object is some- 
what uncertain. Parsamian 21 displays silicate emission 
(Fig. |6]). Integrating the flux of the two components in our 
simple model (Sect. 14. 4p . and correcting the apparent disc 
luminosit y for inclinat i on ef fect {i — 86°, Table [Jjl using 
Equ. 6 of[G reen et al.l l|2006D . we obtained a large relative 
far infrared excess of fd = 75%. These two properties place 
Parsamian 21 into the intermediate-aged category (though 
because of its inclination, i t may actuall y seem younger than 
it is). Following Equ. 7 of iGreen et al. I (|2006. ). from the fa 
value, we also computed the opening angle of the envelope. 
The resulting opening angle of 60° agrees well with the angle 
measured in the direct NACO images (Sect. 1^3. ip . 

Parsamian 21 was placed into the evolutionary scheme 
using two parameters: the silicate feature and the relative far 
infrared excess. In the following we discuss whether its other 
physical characteristics match with those of other FUors. 

(i) Our observations revealed that the circumstellar disc 
of Parsamian 21 is very flat. Due to lack of similar di- 
rect measurements for other FUors, we can only specu- 
late that perhaps all FUors with envelopes have such flat 
discs. On the other hand, the most evolved FUor, FU Ori, 
seems to have no envelope but its disc is probably flared 
jKeiivon Sz Hartmannlll99ll : iGreen et ai1l2006l : iQuanz et"al] 
l2006f ). This might suggest that disc flaring develops at later 
stages, when illumination from the central source may heat 
the disc surface more directly. 

(ii) In a nearly edge-on system like Parsamian 21, one 
expects to see the 10 fj,m silicate feature in absorption. 
The fact that Parsamian 21 has silicate emission indi- 
cates that the line of sight towards the central region is 
not completely obscured. Using the optical depth of the 
15.2 fim CO2 ice f eature, we calcul a ted a n Ay = 8 mag 
(Av = 38.7Ai5.2Mm, ISavage fc Mathid Il979l ). This value 
is surpri singly low compared to V 1057 Cyg {Ay ~ 50 — 
100 mag. lKenvon fc Hartmannlll99lh . This indicates a much 
more tenuous envelope, which is also supported by the low 
envelope mass of 0.02 Mq in our m odelling. 
(Hi) Following lOuanz et al.l |2003), we analysed the profile 
of the 15.2 /xm CO2 ice feature of Parsamian 21. The inset 
in Fig. |6] shows that the feature has a characteristic double- 
peaked sub-struct ure, very similar to HH46IRS, an embed- 
ded young source (|Boogert et aLlbOoj ) . HH 46 IRS is a refer- 
ence case for processed ice. The presence of processed ice in 



Parsamian 21 indicates heating processes and the segrega- 
tion of CO2 and II2O ice, already at this evolutionary stage. 
Other FUors exhibitin g this kind of profil e are LI 551 IRS 5, 
RNO IB and RNO IC ijouanz et al.ll2007l ). 

The evolutionary state of a young stellar obj ect can also 
be est imated following the method proposed bv lChen et al.l 
(|l995l ). According to their Equ. (1) we calculated a bolo- 
metric temperature of Tboi = 410 K for the measured SED. 
We compared this value with the distribution of correspond- 
ing values among young stellar objects in the Taurus and 
pOphiuchus star forming regions (Chen et al. 1995). From 
this check we can conclude that Parsamian 21 s eems to be a 
class I object, and its age is ~ 10^ yr. However. ICreen et al.l 
( 2006) argued that the apparent SED of the disc component 
depends on the inclination. Thus we computed Tboi also for 
a face-on disc configuration and obtained Tboi = 1160 K, 
corresponding to a Class II object. In fact, Parsamian 21 is 
probably close to the Class I / Cla s s II bo rder, in accordance 
with the proposal of lOuanz et all |2003). 



5 SUMMARY 

We present the first high spatial resolution near-infrared di- 
rect and polarimetric observations of Parsamian 21, with 
the VLT/NACO instrument. We complemented these mea- 
surements with archival infrared observations, such as 
HST/WFPC2 imaging, HST/NICMOS polarimetry, Spitzer 
IRAC and MIPS photometry, Spitzer IRS spectroscopy as 
well as ISO photometry. Our main conclusions are the fol- 
lowing: 

(1) We argue that Parsamian 21 is probably an FU Orionis- 
type object; 

(2) Parsamian 21 is not associated with any known rich 
cluster of young stars; 

(3) our measurements reveal a circumstellar envelope, a po- 
lar cavity and an edge-on disc; the disc seems to be geomet- 
rically flat and extends from at least 48 to 360 AU from the 
star; 

(4) the SED is consistent with a simple circumstellar disc 
-I- envelope model; 

(5) within the fram ework of an e voluti onary sequence of 
FUors proposed by lOuanz et all (|2007l ) and |g reen et al.l 
(200d|), Parsamian 21 can be classified as an intermediate- 
aged object. 
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Table Al. Log of IRS calibration measurements. 



Channel 


AOR 


Date 


Target 


Short Low 


16295168 


2005-NOV-21 


HR 7341 




19324160 


2006-Jul-5 


HR 7341 


Short High 


16294912 


2005-NOV-21 


HR 6688 


Long Low 


16463104 


2005-Dec-19 


HR 6606 


Long High 


16101888 


2005-Oct-18 


HR 2491 




Figure Al. Top: Dots show the measured spectrum of HR£ 
obtained with Spitzer/IRS through the Short High slit; the con- 
tinuous line show a model spectrum of the same star (from 
iDecin et al]|2004h . Bottom: Scheme of the data grid; the + sign 
marks the slit centre; dots show the position of the star relative 
to the slit centre; the rectangle indicates the slit itself; the X sign 
marks the position where the spectrum plotted in the top panel 
was taken. 
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APPENDIX A: REFINED CALIBRATION OF 
THE SPITZER/IRS BEAM PROFILES 

In those cases when the source is not well centred in the 
slit, part of the stellar PSF falls outside of the slit, lead- 
ing to flux loss. When the precise absolute flux level of the 
spectrum is important, one should correct for this flux loss. 
Correction for this kind of flux loss is not implemented in 
the Spitzer data reduction pipeline. For point sources flux 
loss in off-centred sources can be efficiently corrected if we 
know the beam profiles of the different IRS slits. In order 
to construct the beam profiles for all four IRS slits, we re- 
duced and analysed dedicated IRS calibration measurements 
taken in spectral mapping mode in a regular grid. Table [XT] 
shows a log of these calibration measurements. We used the 
pipeline processed post-BCD files (pipeline version S14.0.0). 
At each spatial position, we divided the measured spectrum 
with a synthetic MARCS model spectrum appropriate for 
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Figure A2. Beam profiles of the four Spitzer/IRS channels, per- 
pendicular to the slit, at selected wavelengths. Filled dots show 
the measured profile, while open dots indicate the pixelized model 
PSF. 



the measured star (|Decin et al.ll200^ . The result is a data 
cube with two spatial dimensions and one wavelength di- 
mension. We resampled these data cubes to finer spatial 
grid and smoothed them in wavelength. 

We checked whether the measured beam profiles are 
consistent with the PSFs provided by Spitzer's Tiny Tim 
(J. Krist). In Fig. IA2I we plotted the ratio between the ob- 
served and the model spectrum at different distances from 
the slit centre at a certain wavelength. We found slight dif- 
ferences: the measured profiles were in general narrower than 
the model PSFs, and in some cases they were not centred at 
zero. We examined several effects that can cause a difference 
between the model and the measured profiles and we found 
that the differences can be attributed to pixelization effects 
(pixel sizes are 1.8" for SL, 2.3" for SH, 5.1" for LL and 4.5" 
for LH) and to small (few tenths of arcsecond) uncertainties 
in the exact position of the calibration star with respect to 
the slit. Based on our analysis, we decided to use the mea- 
sured profiles to correct our science measurements, and to 
use the Tiny Tim model profiles to estimate the uncertainty 
of the correction. Using this correction, the absolute ffux 
level of IRS spectra has an uncertainty of 10%, but individ- 
ual measurements can be much more precise if the source is 
well-centred in the slit and the necessary correction is small. 
More details on IRS beam profiles will be given in a future 
paper. 

Figure IXSl shows the position of the IRS slits overplot- 
ted on the HST/WFPC2 image of Parsamian 21. While the 
target is well-centred in the Short Low slit, it is slightly 
off-centred in the Long High and it is practically off-slit in 
the Short High. Considering our measured beam profiles, we 
found that no correction is necessary in the Short Low chan- 
nel, a small correction (5 — 20%, depending on wavelength) 
is necessary in the Long High channel, and a significant cor- 
rection must be made in the Short High channel, where the 
corrected fiux is approximately 3-4 times the original one 
(see Fig. IA4|I . After the correction, the different IRS chan- 
nels match fairly well in the overlapping wavelength regimes. 
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Figure A3. Positions of the Spitzer/IRS slits with respect to 

Parsamian 21 (marked with an asterisk). Black rectangle; Short 
Low, small grey rectangle: Short High, big grey rectangle: Long 
High. 




Figure A4. Short High channel spectrum of Parsamian 21. 



